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We consider the Fermi gas in a non-equilibrium state obtained by applying an arbitrary time- 
dependent potential to the Fermi gas in the ground state. We present a general method that gives 
the quantum statistics of any single-particle quantity, such as the charge, total energy or momentum, 
in this non-equilibrium state. We show that the quantum statistics may be found from the solution 
of a matrix Riemann-Hilbert problem. We use the method to study how the finite measuring time 
modifies the distribution of the charge transferred through a biased quantum point contact. 



I. INTRODUCTION 

From a theoretical point of view, a quantum pumpi 
at zero temperature simply excites the Fermi gas from 
the ground state into a non-equilibrium state. The exact 
many particle state obtained after a few pumping cycles 
depends on the applied time-dependent potential. When 
the external potential changes slowly (this is commonly 
referred to as "adiabaticity" in the context of quantum 
pumps) the final non-equilibrium state is completely de- 
termined by the evolution of the scattering matrix at the 
Fermi energy S(t). Therefore, it should be possible to ob- 
tain the distribution functions of various operators (like 
total energy, total momentum, total charge in one of the 
leads, etc) as functionals of S(t). This turns out to be 
a non-trivial problem, which we address in the present 
paper. 

Historically, the first quantity studied in the above set- 
ting was the total energy absorbed by the Fermi gas as 
a result of a sudden introduction of a single s-scatterer. 
The distribution function for absorbed energy shows a 
power-like divergence at low energies - the so called Fermi 
edge singularity (see e.g. Ref.0 for details). The distribu- 
tion of the absorbed energy for a more general situation, 
when the strength of the s-wave scatterer changes arbi- 
trarily in time, was considered In all these cases the 
scattering matrix S(t) commutes with itself at different 
times S(t)S(t') — S(t')S(t). The non-commutative case, 
surprisingly, is much more difficul1i£*£i The first non- 
perturbative treatment of the energy distribution for this 
case& mapped the problem into a one-dimensional chiral 
determinant, and seemed not to be easily adapted for the 
computation of other quantities. In this paper we present 
an alternative approach which is much simpler and also 
readily applicable to distribution functions of any single 
particle operator. 

We illustrate our method using the distribution of 
the charge transmitted through a contact, which has at- 
tracted a lot of attention recently. The second moment 
of the charge distribution (i.e. the shot noise power) for 
a biased quantum point contact was considered theoret- 
ically by Khlus^ and LesovikiS, and the full counting 
statistics was obtained by Lesovik and Levitov (LL) 11 . 



The charge counting statistics for quantum pumps and 
quantum point contact under some specially chosen time- 
dependent bias were also consideredi2iiii4. 

The experimental measurements of the shot 
110 i Ht 4iliiiIiiMi are now well established. The mea- 
surement of the higher cumulants, on the other hand, 
is still an unresolved problem actively discussed in the 
literature^. 

The question of what can really be measured in an ex- 
periment, and the closely related question of the back 
influence of a detector onto the measured system are still 
widely discussed in the Iiteraturai2i22*21. We use concep- 
tually the simplest measurement scheme when the total 
charge Q in the left lead is measured twice. Firstly before 
the pumping is started at time t — 0, and secondly after 
the pumping is finished at t = £/ . At both times the leads 
are decoupled, so Q is a conserved quantity. The first 
measurement projects onto the state with a fixed num- 
ber of electrons in the left lead. After that the pumping 
excites this state into a non-equilibrium state that is no 
longer an eigenstate of the charge operator. The statis- 
tics of the transmitted charge is given by the probability 
to find a particular number of electrons in the left lead 
during the second measurement. 

Alternative approaches discussed in the literature in- 
clude coupling the system to a spin and measuring the 
spin precession angled; or measuring the statistics of the 
photon emission^; or measuring total dissipated energy 6 . 
See also Ref. |^ for a general framework for measuring 
the full counting statistics of a quantum mechanical vari- 
able. 

Our main result is the closed expression for the full 
counting statistics in terms of the scattering matrix 
S(t) and the solution of an auxiliary matrix Riemann- 
Hilbert (RH) problem. We assume that S(i) changes 
adiabatically, but otherwise do not restrict its time de- 
pendence. To demonstrate the power of the method we 
analyze in some detail the case of a biased quantum con- 
tact. In the limit of infinitely long pumping time the 
charge distribution is binomialii. We compute the cor- 
rections to the binomial distribution when the pumping 
time is long but finite and discuss the physical meaning 
of these corrections. 
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The rest of the paper is organized as follows. In sec- 
tion HU we relate the counting statistics to a determinant 
of a certain operator; in section lTTll we compute this deter- 
minant using a solution of a matrix RH problem; in sec- 
tion llVl wc apply the general method to a biased quantum 
point contact and rederive LL's answer. In section Ivl we 
asymptotically solve the relevant RH problem and obtain 
corrections to LL's answer due to finite measuring time. 
Some technical results are relegated to the appendix. 

II. DISTRIBUTION FUNCTION OF 
TRANSMITTED CHARGE 

A. Notations 

The free Hamiltonian Hq describes non-interacting 
Fermi gas in two uncoupled leads. Let the single par- 
ticle states ipne be the eigenvalues of the Hamiltonian Hq 
with energy e, where n = 1 . . . N labels energy degenerate 
states in both the left and the right lead. In the simplest 
case when the leads are one dimensional n = 1 for the 
left lead and n = 2 for the right lead. The annihilation 
operator in the state tp ne is denoted by a ne . We assemble 
N such operators in the TV-dimensional column a e 

/ a u \ 

a e = : (1) 

and use the notation a+ = (a^ e , . . . ,a^ e ) for the TV- 
dimensional row of the creation operators. This gives 
for the free Hamiltonian 

H = ^ea+a £ . 

e 

The time-dependent potential that couples the leads is 
described by an TV x TV matrix M(t,e,e'), so the full 
Hamiltonian has the form 

#(i)=ifo+^a+M(t,e,e'K'- 

Analogously, the total charge operator Q in the left lead 
is given by 

Q = ]Ta+La e , (2) 

e 

where L is the projector on the states in the left lead, 
and charge is measured in units of electron's charge. For 
two one-dimensional leads L = ( J § ) . 

B. The Characteristic Function 

We want to compute the probability P(q) that the 
charge in the left lead changes by q during the time tf. 



The characteristic function x(A) is defined as 

X (A) = ^P(g)e-^. (3) 

At zero bias the many particle ground state |0) corre- 
sponds to filling the states below the Fermi energy in 
both leads 

io>=nrKei>, w 

e<0 n 

where |) is the true vacuum state with no electrons. We 
assume that e takes some discrete values, so the total 
number of electrons is finite, and also set the Fermi en- 
ergy to zero. At zero temperature we have 

X(A) = (0\U + e- tXQ Ue lXQ \0). 

where U = U(tf) is the many particle evolution operator 
which is the solution of the Schrodinger equation 

i^=H(t)U(t), 17(0) = 1 

Since the ground state (QJ is the eigenstate of the oper- 
ator Q, the characteristic function x can be written as 

X = (0|[/+e- lAQ [/|0)(0|e jAQ |0). (5) 

C. One Particle Scattering Problem 

To proceed, we separate the single-particle scattering 
problem from the computation of the expectation value 
over the ground state of the Fermi-gas^. Since the 
Hamiltonian H(t) is quadratic, we have 

L/a+|)=]>>(e, e ')a+|), (6) 

e' 

where <r(e, e') is some N x N matrix that fully describes 
the single particle time-dependent scattering problem. 
We introduce the operator a with matrix elements a(e, e') 
that acts on the full space of all possible states (i.e. in 
both channel space and energy-space). This operator is 
unitary a^ 1 — a + . We show that a can be related to 
the scattering matrix S(t, e) on the instant value of the 
potential M(t) for electron with energy e. 

Consider the scattering states with energy e propagat- 
ing towards the contact ip™ and out of it V'n"*- Let us 
label them and fix their phases in such a way that far from 
the contact i/j ne = ip^ + ipn™* ■> where ip n e is the eigenfunc- 
tion for decoupled electrodes introduced previously. In 
the presence of the potential M(t) the eigen space cor- 
responding to the energy e is still TV-times degenerate. 
Far from the contact an eigen function has asymptotic 
form a n i(>™ + 6 n ^°"*. The scattering matrix S(t, e) 
relates the amplitudes of incoming and outgoing states 
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b n = ^ m S(T,e) nm a m . Note, that the reflection coeffi- 
cients are diagonal matrix elements of S and our choice 
of phases for the scattering states ensures that for decou- 
pled leads 5 = 1. For two one-dimensional leads 



A 

r>* A 

-b 



(7) 



where A is the transmission amplitude and B is the re- 
flection amplitude for the electrons incident from the left. 
Eq implies that the scattering potential is real (i.e. 
there is no magnetic field) and uses this fact to relate 
the scattering amplitudes for the electrons coming from 
the right and from the left. For example, consider the 
two half-infinite one-dimensional conductors occupying 
left and right half-line and separated by an infinite po- 
tential wall of width I. Our choice of scattering states 
dictates A = and B = 1. When the conductors are 
connected and potential wall is completely removed we 
get B = and A = e llp ' n , where p is the momentum of 
electrons with energy e. 

The choice of the scattering states is not unique. An- 
other (more common) approach is to choose the scatter- 
ing states in such a way that the scattering matrix for the 
fully opened channel is one. Denoting the scattering ma- 
trix in this representation by S we see that S = SS~p en 
where S open is the scattering matrix for the fully open 
channel in the representation adopted in this paper. 

The adiabaticity condition (common in describing 
"quantum pumps" ) means that S (r, E) changes slowly 
both as a function of time and energy 



dS" 1 dS 



dt dE 

Using JSJ one can show (see e.g. 

a = S 



< 1. 



[) that 



(8) 



(9) 



when a acts on states near the Fermi energy. Here the 
operator S is the Fourier transform of the scattering ma- 
trix, i.e. it has the matrix elements 



5, 



1 

2tw 



dtS(t)e 



(10) 



where v is the density of states near the Fermi energy 
per single channel. In the rest of the paper we often 
consider operators acting on both energy and channel 
indexes. For example, in addition to a and S introduced 
above we use the operator L that has matrix elements 
L8 ee '. We use hats to distinguish these operators from 
the corresponding matrix- valued functions. The notation 
"Tr" is used for trace in both channel and energy space, 
while the symbol "tr" is reserved to denote the trace in 
the N x N channel space only. 



D. Averaging Over the Fermi Distribution 

The remaining task is to compute the average with 
respect to the Fermi ground state Q}. We combine © 



with J5J) and (0) and introduce the operator 



to obtain 



R = &+e- ixt a = e - lXa+La 



[/+ e -^£/a^|)=5>(e,e')a+|). 



(11) 



The Fermi distribution at zero temperature can be 
viewed as the diagonal operator / with matrix elements 
f(e,e') = 5 ee '9(—e). In the block notation that separates 
the states with positive and negative energies we have 



/ 



1 




In the same block notation the operator R has the form 



R 



Rn R12 

i?21 R22 



(12) 



Since the ground state l@J is simply a Slater determinant 
we have (0|[/ + e~ lAQ L/|0) = det# n for the first term 
in JSJ. This expression is further simplified by noting 
that the operator 1 — / + JR has the same determinant 
as R\\ because in the block notation it takes the form 



Rn R\2 
i 



(13) 



The second term in (0 is simplified because |0) is an 
eigenstate of the charge operator Q 

(0|e lAQ |0) = e ,;A<0|Q|0> = e ,;ATrZ// ". 

In this way we obtain 

X = dct(i-/ + / J R)e 4ATri/ (14) 

Eq. (|TlJl combined with (|TT)) and © gives 

X = det(l - / + /5- 1 e- ai 5)e a ' Iv£/ ". (15) 

In the time representation 5 and L are multiplications 
by the matrix- valued functions S(t) and L respectively 
and / has matrix elements 



/(M') = 



1 



2tt t - V + iO 



(16) 



The representation of counting statistics for a biased 
quantum point contact through the determinant of type 
was obtained in Ref. and later adopted for para- 
metric pumpingl 3 and re-derived in different ways^. Our 
derivation is close in spirit to Ref>i£ and is probably the 
simplest. In addition it leads to a simple regularisation 
procedure which we discuss below. 
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E. Biased Quantum Pump 

Finally, we briefly mention that in the case when the 
left lead is biased with respect to the right lead by a time- 
dependent potential V(t), a gauge transformation can be 
performed on all states in the left electrode 

a fc — > e Jo y ' ak- 
in the new gauge the ground state is again given by 10} . 
The distribution function for this case can therefore be re- 
covered from HI 5(1 by substituting the gauge transformed 
scattering matrix 

e iLftV(T)dT Se -zLfZV(r)dT ( 17 ) 

instead of S. 



The first term in Eq 1(20(1 contains the information about 
the second and higher moments of the charge distribu- 
tion. 

III. RELATION TO A RIEMANN-HILBERT 
PROBLEM 

In this section we relate the operator trace 

Tr = Tr (ln(l - / + f R) - f In Ft) (22) 

from Eq. 1(20(1 to a solution of an auxiliary N x N matrix 
Riemann-Hilbert problem. 

A. General case 



F. Regularisation 



We use 



Eq. © is only valid for the states near the Fermi 
surface. Since the determinant in 1 (14(1 contains a non- 
vanishing contribution from the states deep below the 
Fermi surface, the transition from 1(14(1 to ((15(1 is not 
quite correct. In particular, the determinant in ((15(1 is 
formally infinite (and so requires regularisation) while 
Eq. 1(14(1 gives a finite answer. To find out the correct 
regulator we use the decomposition 

ln X = Tr(ln(i-/ + /i?.)-/lni?) (18) 

+ Tr(f]nR) + iXTrLf. 

In the first term in ((18(1 both the contributions from the 
states deep below the Fermi surface (when / = 1) and 
the states high above the Fermi surface (where / = 0) are 
zero, so Eq. (JHJ can be used. The two last terms in 118(1 

can be re-arranged to give — iX Tr facr, f]a + L) = —iXQ av 

where the square brackets denote the commutator. Again 
the contribution from the states away from the Fermi 
energy vanishes and we obtain 

Qav = Tr f [S, nS-'L) = ±. [ tri^-S- 1 ^ (19) 



2tt 



dt 



To obtain the last equality in ((191) we compute the trace 
in the time representation, where all operators, except / 
are diagonal, and matrix elements of / are given by ((16(1 . 
Comparing ((19(1 with Brouwer's formula™ for the av- 
erage charge transmitted through a quantum pump, we 
identify Q av with the average value of the total transmit- 
ted charge. 

Combining these results together we get 



ln% = Tr 



w 



1-f + fR 



iXQ av . (20) 



where in the time representation R is multiplication by 
the matrix valued function 



R(t) = S- 1 (t)e- aL S(t). 



(21) 



Tr = 



dXTr 



CI-/ + /V?! 7- ■ )^ 



(23) 

To invert 1 — / + fR we need to solve the following aux- 
iliary RH problem. Find a matrix-valued function Y(t) 
of the complex variable t such that Y 

• is analytic on the complement of the cut (0,£/) 

• obeys 

Y-(t)y+ 1 (t) = R{t)e lXL when t G (0,t/), (24) 
where Y±{t) = Y(t±i0) 

• tends to 1 at infinity 

Y — ► 1 when \t\ — > oo. 

Using 1(16(1 and analyticity of Y± we observe that Y± has 
the following properties 

fY+f = Y+f. 

Using the above and substituting R = YJY+ x e~ lXL 
one can verify 

(l - / + fR) ^ = e iXL Y + ((1 - f)Y^e-* XL + fYl 
which after substitution in 1(231) gives 



Tr 



dX 
2^ 



dt tr e 



Combining ((25(1 with ((20(1 and 119(1 we obtain our main 
result 



f x dX 
lnX = / 27 



dt 



(26) 



tr ■ 



' d (Se iXL Y + ) 



dt 



(S, 



iXL 



Y+ 



L 
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where Y is analytic in the complement of the cut (0, tj), 
obeys = S~ 1 e~ iXL Se iXL along the cut and goes 

to 1 at infinity. 



B. Charge Transfer per Cycle for Periodic 
Pumping 



where 



|A|V A + \B\ 2 



P = 



A*B(1 



(32) 



e lX ) 



-AB*(l-e 



Eq. I|26|) can also be used when the scattering matrix 
is periodic S(t) = S(t + T) and we are interested in the 
limit of the large number of cycles (i.e. tt/T 3> 1). The 
RH problem in this case can be somewhat simplified by 
rolling the infinite real axis into a unit circle with the 
mapping z = e 27Tlt / T . After that the RH problem J2U 
reduces to finding two functions Y± (z) such that 

• Y + is analytic for \z\ < 1 and Y_ is analytic for 
\z\>l 



• on the unit circle \z\ = 1 



y_y: 



1 when 



S- 1 (z)e- lXL S(z)e 



i\L 



(27) 



The distribution of charge transmitted during time tj 
in the limit tj/T 1 is now given by 



(28) 




i\L 



Y. 



L 



IV. COUNTING STATISTICS FOR A BIASED 
QUANTUM POINT CONTACT 

In this section we apply the technique introduced above 
to the quantum point contact described by the scattering 
matrix (JJJ and biased by the time-independent potential 
V. 



A. General formulas 

After the gauge transformation (|17|l the scattering ma- 
trix takes the form 



S(t) 



JLVt 



S Q e 



-iLVt 



(29) 



where the time- independent matrix (7J is denoted by So- 
The matrix S~ 1 e~ lXL Se lXL can be represented as 



s -l e -i\L Se i\L 



e lLVt R e- lLVt (30) 



lXL Soe lXL is time- independent. It is 



where i?o = Sq 1 ( 
convenient to decompose Rq into a product of upper and 
lower-triangular matrices 



ft? 



1 (3 
1 



(31) 



B. Long time limit i/F> 1 

The scattering matrix l|29|) is periodic with period T = 
2?t/V and in the limit tfV > 1 we can apply the results 
of section IfflBl The RH problem (gjjl takes the form 



Y-Y7 1 = 



1 

« 1 



a 
h 



1 pz 
1 



(33) 



where z 



and we used l|3(J|l and Q3ip. 



The solution of RH problem (|33() is obvious now 



Y- 
Y, 



1 
— 1 

Z 

1 




-pz 
1 



a 



Substituting Y + into 128|) we get 



lnx(A) 



tjV 
2ir 



In a 



(34) 



with a given by (|32|l . The above formula is LL's answer 11 
for the counting statistics of a biased quantum point con- 
tact. Thus we have verified that measuring charge in 
one lead after the leads are disconnected yields the same 
counting statistics (in the limit tfV 1) as the measure- 
ment scheme based on a coupling to a spin, as suggested 
by LL. 



V. CORRECTIONS TO THE COUNTING 
STATISTICS DUE TO FINITE MEASURING 
TIME 

A. general considerations 

Eq (|34|l gives the leading term in the expansion of In x 
in tfV. The correction is proportional to IritfV, and to 
compute it we need to solve RH problem (1241) 



Y-YZ = 



1 

,-ivt 1 



a 
± 



1 Pe lVt 
1 



(35) 



This RH problem is no longer periodic, since we need 
to take into account explicitly the finite length of the 
cut. We are not aware of an analytical solution for RH 
problem (|35|l . Instead, we find an approximate solution 
in the limit tfV 3> 1. To do this we use the saddle- 
point method introduced for RH problems in Ref. |U 
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The idea is to modify RH problem l|35|) in such a way 
that the cuts become vertical and the oscillating factors 
e ±iVt t urn jp^Q decaying ones. To implement the plan, 
we construct an auxiliary matrix- valued function W and 
look for the solution in the form Y = W£. It turns out 
that £ obeys another RH problem, where all cuts are 
vertical, and the jump functions decay along the cuts as 
e~l*l v '. As a result £ is close to 1 away from the branch 
points and tf, and W yields the desired approximate 
solution to RH problem (J3SJ. 



© 



FIG. 1: Regions of analyticity for W[t). 



B. approximate solution of RH problem (1351) 



It is convenient to introduce the function 
In a 



ip(z) = exp 



2„ 7f<7 



(36) 



which has the same region of analyticity as Y, tend to 1 
at infinity and obeys 



a 
i 



(37) 



along the cut (0,i/). 

The approximate solution W(t) is constructed as ana- 
lytic in the complement of three cuts. The first two cuts 
are along the two vertical lines Re t = and Re t = t / and 
the third cut is along the real axis between and t / (see 
fig ^| . The cuts divide the complex plane into four re- 
gions : the half-plane Ret < (region 1), the upper part 
of the strip strip < Ret < tf (region 2), the lower part 
of the same strip (region 3) and the half-plane Re t > tf 
(region 4) . Denoting the value of the function W inside 
the region i as Wi{t) we define W in the following way: 



Wi = w 3 



vt 



1 (3d 
1 

1 n 



Comparing l|37|) and (|31|l with l|35|l we note the important 
property 



W-W+ 1 = Y-Y+ 1 



(38) 



along the cut (0, tf). 

We are now looking for the solution of the original RH 
problem l|35|) in the form Y = W£. Combining (|35|l with 
the definition of W we conclude that £ is analytic in each 
of the three regions Ret < (i.e. region 1),0 < Ret < tf 
(this is regions 2 and 3 combined together) and Ret > 
tf (i.e. region 4). Due to condition (|3"5jl the function £ 
does not have a cut between regions 2 and 3. The other 
conditions on the function £ come from the requirement 
that Y is continuous on the vertical cuts of W. Using 



the notation £j for the value of £ in the region i, this 
condition gives on the upper half of the cut Re t = 



= W^W 2 



1 (3e lVt 
1 



(39) 



and on the lower half of this cut 



ac 1 = wr'ws = ( ae iv* f) 4>. (40) 

Importantly, the right hand sides of Eqs. (|39|l . (|40|) are 
exponentially close to 1 when t goes to infinity along the 
cuts. The difference is of the order e~'*' y and can be 
neglected when \t\ ^s> 1/V. The same happens along the 
other cut at Re t = tf. Using this fact alone we prove (see 
appendix that the function is close to 1 when t is 
far away from and t * . More precisely, in the region 



\t\ > and \t 



(41) 



£ is close to 1, so Y is close to W, as required. Explicitly, 
we obtain for Y + 



Y + (t) 



-Pe lVt 
1 



when 
when 
when 



t < 

< t < t f 

t f <t 



(42) 

where ip is given by l|36() . Eq. I|42() is not valid in the 
vicinity of the end points and t * , where conditions l|41|) 
are not met. It turns out that we do not need the exact 
behavior of Y+ in those regions to obtain the hitfV term 
in ln%. 



Computation of the distribution function and 
interpretation of the result 



The function Y + (t) enters ln^ through the combina- 
tion ^1"+ 1 (see Eq. 123 • There are, therefore, two 
contribution to ln^- The first contribution comes from 
differentiating the factor e lVt in (|42|l and reproduces 
Eq. (|34|l which is linear in tfV. The second contribu- 
tion comes from differentiating ip in (|42|l and is propor- 
tional to IntfV. Combining both contributions together 
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we obtain 



lnx = 



tjV 
2vr 



111, 



hxtfV 
2tt 2 



In 2 a, 



(43) 



where a is given by i|32|l . 

Finally, there are contributions from the regions t < 
l/V and \t — tf \ < l/V where the asymptotic solution 
constructed above is not valid. These contributions are 
of order ln(rV), where r is the switching time. Assum- 
ing t ^ 1/V these extra logarithms can be neglected. 
Combining all the conditions specified above we find the 
validity region for Eq. I|43|) 



l/V « t f . 



(44) 



The characteristic function l|43|) corresponds (with log- 
arithmic accuracy) to the probability distribution 



P{k) =5^p(n)B(n,fc), 



(45) 



which is a convolution of a Gaussian distribution for the 
number of attempts 



p(n) oc e 



V(n-^) 2 /(41n^) 



(46) 



and the binomial distribution 



B(n,k) 



k\(n - k)\ 



\A\ 2k \B\ 2i - n ~ k \ (47) 



for the number of successes. The mean number of at- 
tempts is given by tf V/2tt and the variance is logarithmic 

2 l tjV 

er ~ In -4r~ ■ 

There are, therefore, two sources of shot noise in a 
quantum point contact. One source is the fluctuations 
in the number of incident electrons (number of attempts 
to go through the contact) and the other is the fluctu- 
ations in the number of reflections (number of failures). 
According to (|4*5l) these two sources are statistically in- 
dependent. The physical picture of the two statistically 
independent sources of quantum shot noise was suggested 
(but not proven) by LLAi. Note, that in order to confirm 
the statistical independence we needed to obtain the sec- 
ond term in Eq. I|43[) and analyze its dependence on the 
transmission amplitude A. 



VI. DISTRIBUTION FUNCTION OF THE 
DISSIPATED ENERGY 

We mention without derivation that the characteristic 
function Xe(A) = / P(E)e- lXE dE for the distribution of 
energy dissipated during the pumping^, is given by the 
expression 



InXe(A) = — J d\ 
— iXE av , 



(48) 



which is very similar in form to Eq. I|26[l . The analog of 
the Ricman-Hilbert problem (|24ll is Y- e Y+} = Re, where 



R e (r) =S- 1 (t)S(t + X). 
The average dissipated energy 2 ^ 



1 f tf JSdS' 1 



E av = ^l tr( 



' dr dr 



)dr. 



(49) 



(50) 



appears analogously to Q av during the regularisation of 
the determinant. 

Similarly, the distribution function of the values of any 
single-particle operator is given by expression similar to 
Eq. J2EJ), with appropriate replacement of R and 



VII. CONCLUSION 

It is clear from the derivation, that the charge operator 
Q in equation Q can be replaced with any single-particle 
operator. This will modify the counter-term l|19|) . as well 
as the relation (|21|l between R and S. Also Eq. I|26l that 
relates the characteristic function with the solution of the 
RH problem is modified. On the other hand, the general 
form of the answer remains the same, i.e. the distribution 
function is still expressed as an integral over an expres- 
sion containing the solution of a particular matrix RH 
problem. 

The general method is used to compute the distribu- 
tion of the charge transmitted through a biased quantum 
point contact in finite time. We confirm LL's picture of 
the two sources of the charge fluctuation (see section lV C|) 
and prove that the transmission attempts are statistically 
independent from the reflections when the total observa- 
tion time is much longer than inverse bias. 



APPENDIX A: ASYMPTOTIC EXPANSION 
FOR £ 



To complete the discussion in section IV Bl we need to 
find the function £(t) with the following properties : 

• £ is analytic in the complement of the two vertical 
cuts Re t — and Re t = tf 



• £ obeys 



= w^Wt 



(Al) 



along the cut Ret = 0, where i — 2,3 and notations 
are the same as in section IV Bl 



• £ obeys £4^ 1 = W 4 1 Wi along the cut Ret = tf. 

• £ — > 1 when \t\ 



00 



8 



In this appendix we show that £ is close to 1 in the limit 
Vtf 3> 1 everywhere except in the two small regions near 
t = and t = tf 

The general method can be illustrated by finding a 
function r\ with just one vertical cut along the upper part 
of imaginary axis (Ret — Q,lmt > 0) obeying the con- 
dition r\\r\2 X = W^ x Wi along that cut. The case with 
two cuts is exactly the same. 

Analyticity of r\ together with the boundary condition 
at infinity leads to the representation 



r\(z) 



1 

2~7ri 



r)(t'- 0) - <n(t' + 0) 



Rct'=0;lmt'>0 



t' -z 



(It. 



Using the above we obtain the integral equation for the 
values of 771 (t) on the line Re t = 



m{t) = 1 + 

j_ r 

2ni 



(A2) 



'Rct'=0;lmt'>0 



t'-t + Q 



Because of exponential decay of the kernel along the cut 



the integral in Eq (|A2I) converges at t' ~ 1/V -C £, so in 
the limit > 1 we can neglect t' in the denominator 
and get the expansion 



where 



1 C 

m = 1 H 

1 Vt 



C 



V 



'Ret'=0;lmt'>0 

is some constant. 



(i-w 2 - 1 (0^i(* , ))m(t'M* / ~i 
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